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Abstract 

The process e + e~ — > ttH is considered in the kinematic end point region where the Higgs energy 
is close to its maximal energy. In perturbative QCD, using the loop expansion, the amplitudes are 
plagued by Coulomb singularities that need to be resummed. We show that the QCD dynamics 
in this end point region is governed by nonrelativistic heavy quarkonium dynamics, and we use a 
nonrelativistic effective theory to compute the Higgs energy distribution at leading and next-to- 
leading-logarithmic approximation in the nonrelativistic expansion. Updated umbers for the total 
cross section including the summations in the Higgs energy end point region are presented. 
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I. INTRODUCTION 



It is one of the most important tasks of future collider experiments to unravel in detail the 
mechanism of mass generation and electroweak symmetry breaking. Within the Standard 
Model of particle physics (SM) electroweak symmetry breaking is achieved by the Higgs 
mechanism which postulates the existence of an electrically neutral scalar field that interacts 
with all SM particles carrying non-zero hypercharge and weak isospin. The particle masses 
are then generated by the Higgs field vacuum expectation value V = (V2G F ) 1/2 « 246 GeV, 
Gf being the Fermi constant, that arises through the Higgs self interactions. The mechanism 
also predicts that the Higgs field can be manifest as a massive, elementary, scalar Bose 
particle that can be produced in collider experiments. The mass of the Higgs boson is 
expected to lie between the current lower experimental limit of 114.4 GeV |l[ and about 
1 TeV. Current analyses of electroweak precision observables yield a 95% CL upper bound of 
260 GeV for the Higgs boson mass • While a Higgs boson with a mass smaller than 1 TeV 
can be found at the LHC, precise and model-independent measurements of its quantum 
numbers and couplings can be gained from the e + e~ Linear Collider 0,0,0]. 

An important prediction of the Higgs mechanism is that its Yukawa couplings to quarks 
X q are related to the quark masses by m q = X q V. This makes the measurement of the 
Yukawa coupling to the top quark particularly important since it can be expected that it 
will be determined with the highest precision among the Yukawa couplings, either through 
direct or indirect methods. At the e + e~ Linear Collider the top quark Yukawa coupling can 
be measured from top quark pair production associated with a Higgs boson, e + e~ — > tiH. 
This process is particulary suited for a light Higgs boson since the cross section is then 
large and the rate is dominated by the amplitudes describing Higgs radiation off the top 
or the antitop quark. Assuming an experimental precision at the percent level, QCD and 
electroweak radiative corrections need to be accounted for in the theoretical predictions. 

Approximating the top quark and the Higgs boson as stable particles the Born cross 
section was already determined some time ago in Refs. @]. For the 0{a s ) QCD one-loop 
corrections a number of references in various approximations exist. An analysis for large 
energy and small Higgs masses can be found in Ref. pj, while a full computation of the 
numerically dominant virtual photon exchange contributions, where the Higgs is radiated 
exclusively from the top and antitop quark, was presented in Ref. Finally, the full 
set of QCD corrections was given in Ref. 0. On the other hand, the full set of one-loop 



electroweak corrections was obtained in Refs. jlOL and also in Refs. [12]. In Ref. [TJ| a 
detailed analysis of various differential distributions of the cross section a(e + e~ — ► tiH) was 
given. 

A particularly interesting kinematical phase space region is where the energy of the Higgs 
boson is large and close to its kinematic end point. In this end point region the ti pair is 
forced to become collinear and to fly opposite to the Higgs direction in order to maximize the 
momentum necessary to balance for the large Higgs momentum, see Fig. \I]for an illustration. 
The general relation between the Higgs energy Eh and the ti invariant mass Q 2 = (k± + k 2 ) 2 
reads 

Eh = (s + m 2 H -Q 2 ) , (1) 

where \fs is the cm. energy and m# the Higgs boson mass. Thus for large Eh the ti 
invariant mass approaches 4m 2 and the top quark pair is nonrelativistic in its center-of- 
mass system. Since a light Higgs boson having a mass below the W + W~ threshold is quite 
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FIG. 1: Typical constellation of momenta for the process e + e~ — > ttH in the large Higgs energy 
end point region. 



narrow 1 and the hadronic Higgs decay final state factorizes to a very good approximation, the 
strong interactions between the ti pair and the hadronic Higgs final state can be neglected. 
Therefore, close to the Higgs energy end point, the ti QCD dynamics is exclusively governed 
by the nonrelativistic physics known from the process e + e~ — > ti in the ti threshold region 
at y/s » 2m t Q. 

In this regime the so-called Coulomb singularities cx (oc 3 /v ) n , with v — (1 — 4m 2 /Q 2 ) 1 / 2 
being the top quark relative velocity in the ti cm. frame, arise in the amplitudes and render 
the perturbative expansion in the number of loops inapplicable. This singularity structure 
is most easily visible in the Higgs energy distribution, da(e + e~ — > tiH) / dEn- While the 
Born distribution approaches zero for Eh — » -E# ax , da/dEn ~ v [aL the 0(a s ) fixed- 
order perturbative corrections are proportional to a s at the end point [allj and the 0{a 2 s ) 
corrections even diverge like a 2 s /v . In principle, the problem can be avoided by imposing a 
cut on the Higgs energy (or the it invariant mass Q 2 ), but such a measure is unnecessary 
because there exists an elaborate technology from the threshold region in the process e + e~ — > 
ti that allows to sum the Coulomb singularities to all orders in a s and to carry out a 
simultaneous expansion in a s and v using a nonrelativistic effective theory. Imposing a cut 
is also disadvantageous since the nonrelativistic portion of the tiH phase space increases 
for smaller cm. energies or larger Higgs masses. Due to the large top quark SM width 
r t w 1.5 GeV, which serves as an infrared cutoff, the corresponding QCD effective theory 
computations can be carried out with perturbative methods for all Higgs energies in the end 
point region. The effective theory also allows for a systematic summation of logarithmic 
terms oc (cv s In v ) n to all orders in a s using the velocity renormalization group [loT ]^ 

In this paper we compute the Higgs energy distribution da/dEn in the large Higgs energy 
end point region at leading-logarithmic (LL) and next-to-leading-logarithmic (NLL) order in 
the nonrelativistic expansion in the framework of "velocity" nonrelativistic QCD (vNRQCD) 
using the conventions and notations of Refs. (TH [3, 17, [3, [HI]. 2 We neglect NLL order 



1 For m H = 115(150) GcV one finds T H = 0.003(0.017) GeV 0. 

2 For an alternative approach, see Refs. [20j. 
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effects coming from the top quark finite lifetime that were discussed recently in Ref. |21 
Our analysis contributes to an improved understanding of uncertainties of higher order 
QCD corrections in the large Higgs energy region and is phenomenologically important for 
the total cross section if the e + e~ cm. energy is not too large. 

The program of the paper is as follows: In Sec. |H] the ingredients of the effective theory 
that are necessary for the calculation at hand are presented. Within the effective theory 
framework the LL Higgs energy distribution in the large Higgs energy end point is calculated 
explicitly as an illustration in Sec. lIIIl The computation of the NLL Higgs energy distribution 
including the numerical determination of the NLL matching conditions is outlined in Sec. II VI 
A numerical analysis of the results is performed in Section and Section IVII contains the 
conclusion. 



II. EFFECTIVE THEORY INGREDIENTS 



In our calculations we employ vNRQCD ("velocity" NRQCD), an effective theory which 
describes the nonrelativistic dynamics of top quark pairs where the nonrelativistic scales 
m t v (momentum) and m t v 2 (energy) are larger than the hadronization scale Aqcd- It has 
become the standard convention to call the momentum scale "soft" and the energy scale 
"ultrasoft". In the following we summarize the ingredients necessary for the description of 
the nonrelativistic ti dynamics at the NLL order approximation in the ti cm. frame. For 
details on the conceptual aspects, concerning powercounting, the operator structure of the 
effective theory action, and renormalization we refer to Refs. H EE El El 13 ■ 

The particle-antiparticle propagation is described by the terms in the effective theory 
Lagrangian that are bilinear in the top quark and antitop quark fields, 



— 1 — — 5m t 

2m t 2 



(2) 



where the fields ip p and Xp destroy top and antitop quarks with soft three-momentum p in 
the ti cm. frame and T t is the on-shell top quark decay width. The x-dependence of the 
fields describes ultrasoft (~ m t v 2 ) fluctuations while the labels p refer to the soft (~ m t v ) 
three-momentum of the quarks. The term 5m t is a residual mass term that is specific to the 
top quark mass definition that is being used. In order to avoid the pole mass renormalon 



problem (22I . 23l | we employ the IS mass scheme 0, |25| with 



5m t 
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(5) 
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where j3o = 11/3(7^— 4/3Tn/ is the one-loop QCD beta function, ai = 31/9Ca— 20/9Tn/ the 
coefficient of the one-loop correction to the effective Coulomb potential and Ca = 3, Cp = 
4/3, T = 1/2 are SU(3) group theoretical factors. For the number of light flavors we take 
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rif = 5. The parameter v is the dimension-zero vNRQCD renormalization scaling parameter 
used to describe the correlated running of soft and ultrasoft effects in the effective theory 
governed by the velocity renormalization group jl^]. 3 The resulting top/antitop propagator 
in the effective theory has the form 



pO _ pi J (2m t ) + iY t /2 — 5m t 



(7) 



Up to NLL order the top-antitop quark pair interacts only through the effective Coulomb 
potential 00, 



V c (p,q) 



V C (S) M 



4nC F a s (m t v) \ a s {m t v 



k A k' 

V<f\v) = -AnC F a s (m t u) , 



47T 



-A) In 



(8) 



where k = p — q is the momentum transfer. The term Vc S \v) is the color singlet Wilson 
coefficient of the 4-quark Coulomb potential operator. Here, v = 1 corresponds to the 
hard scale at which the effective theory is matched to the full theory, and v — Vq, Vq being 
of the order of the typical tt relative velocity, is the scale where the matrix elements are 
computed. The evolution of the Wilson coefficients from the matching scale down to the 
low-energy scale sums logarithms of the velocity to all orders and is governed by the velocity 
renormalization group equations |l5[ which are determined from the anomalous dimensions 
of the effective theory. 

Top-antitop quark production in the nonrelativistic regime in the LL and NLL approxi- 
mation in a 3 5*i spin triplet or a 1 5'o spin singlet state is described by the currents 



J l P = ^l^ji^X-p 



(9) 



where Cij(u) and Cq{v) are the corresponding Wilson coefficients. The currents do not 
run at LL order, but they have a non-trivial anomalous dimension at NLL order from UV 
divergences in effective theory two-loop vertex diagrams 0> E 0] • The NLL running of 
the Wilson coefficients reads 



cij(f) = cij(l) exp(/(z/,2)) 
cbH = cq(1) exp (f(y,0)) , 



(10) 



The renormalization scales for soft and ultrasoft fluctuations, /is and fijj, are correlated through the heavy 
quark equation of motion, fijj = u%/mt. This correlation is necessary for a consistent renormalization 
of the nonrelativistic effective theory [26|. The correlated running from the hard scale down to the soft 
and ultrasoft scales is described by the dimensionless scaling parameter v defined by [is = rntv and 
fMjj = mtv 2 . Thus, v = 1 corresponds to the hard matching scale, and effective theory matrix elements 
are evaluated for v w 0.2, i.e. of order of the typical top quark velocity. 
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where 



f{v,S 2 ) = +a 2 -na s {m t ) (1 - z) + a 3 na s (m t )ln(z) 



+ CZ4 ira s (m t ) 



1 _ z l-13C A /(6po) 



+ a 5 7!-a s (m t ) 



1 _ z l-2C A //3o 



+ a vra; s (mt)) (2 — 1) — w 1 ln(w) 



(11) 



with 



and 



fl 2 



«3 



Q s (m f v) 
a s (m t ) 



w 



a s (m t v) 



C F [C A C f (9C a - IOOCf) - (3 (26C 2 A + 19C A C F - 32C 2 F )] 



26 (3 2 C A 



(14 



a 



(3 2 (Q(3 -13C A )(P -2C A 
-Q(3 2 \2C F + C A (S 2 - 3) 

24C 2 (11C A -3(3 )(5C A +8C F ) 
13C A (6(3 -13C A ) 2 

8C F (C A + C F ) (C A + 2C F ) 
3f3 2 



C 2 A (9C A - 100C F ) + Po C A 74C F + C A (13S 2 - 42) 



a 5 



C 2 [C A (l5-US 2 )+f3 (4S 2 -3 
Q((3 -2C A ) 2 



(12) 
(13) 



The term S 2 is the square of the total ti spin. We take the convention where the matching 
conditions at v — 1 only account for QCD effects, so at LL order we have ci(l) = Co(l) = 1. 

Through the optical theorem the ti production rate for a ti invariant mass Q 2 ~ Am 2 
involves the imaginary part of the time-ordered product of the production and annihilation 
currents defined in Eqs. Q, 

<(Q 2 ,m 4 ,z/) = iJ2 j d A xe- l ^ x (^0\t jl pl {0)Jt P {x)\o) = 2N c 6 lk G c (a, v, m t , v) , (14) 



p,p 



x e 



A (Q 2 ,m t ,v) = Jd 
p,p' 

where q = (y^Q 2 — 2m t , 0) and 



-iq-x 



TJl pf (0)J , P (x) 



N c G c (a,v,m t ,u), (15) 



Q 2 — 2m t — 25m t + iT t 
m t 



(16) 



is the cm. top quark relative velocity. The term G c is the zero-distance S-wave Coulomb 
Green function of the nonrelativistic Schrodinger equation with the potential in Eq. (jHJ. At 
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LL order (i.e. including only the first term on the RHS of Eq. (Qfl) ) the Green function has 
a simple analytic form and reads (in dimensional regularization) [19| 



Gl L (a,v,m t ,v) = — 



liv — a 


In ( 













1 f ia 
- + ln2 + 7 £ + ^ l 

2 ' ^ V 2v 



mf a 1 

H 

47r 4e 



(17) 



Explicit analytic expressions exist for the NLL order correction to G c coming from one 
insertion of the 0(a 2 ) contributions to the effective Coulomb potential in Eq. (jHJ) |3Qj . We 
use the numerical techniques and codes of the TOPPIC program developed in Ref. j3jj 
(see also Ref. [32|) and determine an exact solution of the full NLL Schrodinger equation 
following the approach of Refs. [l^ . 



III. HIGGS ENERGY DISTRIBUTION AT LL ORDER 



To illustrate the method of computing the Higgs energy distribution for large Higgs 
energies we start by considering the simplified case where the process is mediated through 
a virtual photon only. The amplitude for the process 



e+(k')e-{k) - 7 * -> t{h)i{k 2 )H{ PH ) 

in the full theory reads 

AT\Q e Q t a 



(18) 



M~f = i\ t 



' v e +{k)^ ll u e -(k')} 
i>H + kfi+ m t 



1 \(ph + h) 2 - m 2 t 



t + r 



-j>H - ^2 + m t 
(Ph + h) 2 ~ m 2 



Vt(k 2 ) 



where 



e m t 
2s w M w 



(19) 



(20) 



e being the electric charge and s w (c w ) the sine (cosine) of the Weinberg angle. In the large 
Higgs energy region one has 



E 



H 



1 



kt « k 2 « kt = (E t ,kt) 



(s + m 2 H — Am 2 ) 



m? + IkJ 2 



t\ j 



1 / ( g + m^-4m 2 ) 2 2 
1 11 2 I 4s 



1/2 



where 



/\2 



(fc + fc') 



(21) 



(22) 
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and one can relate the effective theory top and antitop spinors in the tt cm. frame with the 
full theory spinors by a Lorentz boost, 



ut{h) = U(k t 







E t + m t 
2m t 



1/2 



Et+mt 



Z7 



U(k t 



X—p 




E t + m t 
2m ( 



1/2 



Et+mt 



{W2XI 



(23) 



where p is the small residual top quark momentum in the ti cm. frame. Note that on the 
RHS's shown in Eqs. (}23|) we have only kept the effective theory operators that are leading 
order in the expansion for p <C m t in the ti rest frame. This is sufficient for the matching 
computation at LL and NLL order. We also note that in the effective theory we describe 
the antitop quark by a positive energy spinor which makes the charge conjugation operation 
in the second equation necessary. In the large Higgs energy end point region the amplitude 
then reduces to the form 



AnQ e Q t a 
s 

1^~s 



v e +(k')^iU e -(k) ^lcr j (ia 2 )x- p 
E t 



rplj 



(Ph + h 



mi 



m t 



1.1 



m t {E t + m t ) 



(24) 



We see that for virtual photon exchange the top-antitop pair is produced only in the spin 
triplet configuration. Without QCD corrections, the distribution in the large Higgs energy 
region then reads 



da 



2s(2vr) 5 4 ^ 

v spins 



Ph, 



{2mt? lE \M J {E H *E%r)\ 2 6n(k + k>-k 1 -k 2 -. 

x 5+(p 2 H — m 2 H ) 5+(k\ — m 2 ) 5+(k 2 — m\) d 4 pn d 4 ki d 4 k 2 



57TS 



5/2 



Lij Tik T* t Im 



x 



A[ k ((k + k' -p H )\m t ) 
[(1 + x H - Ax t f - Ax H ] 1/2 dE H dVL H 



2s 5' 1 - 



k { k j 



(25) 



where 5+ indicates that in the on-shell condition only positive energies are accounted for. 
The term A lk is the current correlator in Eq. (jl4)l without QCD effects and for stable top 
quarks, Ly the lepton tensor and 



m. 



x t 



■''■I = 



m 2 H 



mi 



x z = 



(26) 



In the large Higgs energy region, including the QCD effects coming from the Coulomb 
potential in Eq. (jHJ) and the finite top quark lifetime, the LL result for the Higgs energy 
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distribution reads 



8 N F^ 1/2 
da = 3/2 c ^ 2 Im [ Gl L (a, v, m t: v) ] [(1 + x H - 4x 4 ) 2 - Ax H ] dE H 



\ 2 QlQ 2 a 2 (l-x H + 4x t ) 2 + 8x t 
3 (l + x H - 4x t y 



(27) 



If the virtual Z-exchange is accounted for, the computation is slightly more involved since we 
also have to account for the top pair production in the singlet spin state configuration and 
both correlators in Eqs. ()14I15|) are needed. Including the Wilson coefficient of the currents 
in Eqs. (J10J) as well, the result can be written in the form 

da -{E H *E^) = 8iYc[(1 + ^: 2 4 f 2 - 4 ^ ]1/2 U^F^c\iy)Fr 



x lm[G c (a,v,m t ,u)} , (28) 



where 



p7 ,z _ « 2A ? + 4x t ) 2 + 8x t f 2 2 v 2 (a? e + vp 2Q e Q t v e v t 



l + x//-4a; i ) 2 V 6 (l-x z ) 2 1 - x z 

Z 1 (1 + Xtf - Ax t )(Ax t - x z )(l - x z ) \ (l-x z ) tet 



a 2 g 2 Z v 2 t (a 2 e + gg) ((1 - + 4x f ) 2 + 32s f )a: z 
24 (4x<-xz) 2 (l-xz) 2 ' 



a 2 g 2 z a 2 (a 2 + w 2 ) (1 - x H + 4x 4 ) 2 - 16x t 
"24 (l-x z ) 2 x z 



= ^ 71 ' W 



and 



_ T/ - 2Q /g 2 T( e 

^°w^w ^°w^w ^ W^W 

Here Qf is the charge for fermion /, and T[ is the third component of weak isospin. In 
Fig. it is shown that the singlet contribution governed by the coefficient Fq is negligible 
for cm. energies around 500 GeV. For larger cm. energies it becomes more important, but 
remains smaller than the triplet contribution for all available energies. Figure El also illus- 
trates that the triplet contribution, determined by the coefficient F^' , is dominated by the 
term proportional to the square of the Yukawa coupling X t . This term describes the effects 
of Higgs radiation off one of the top quarks. 

We note that the NLL (0(a s )) matching conditions for the three triplet Wilson coeffi- 
cients Cij depend on the ti spin configuration (i.e. on j) since the kinematic situation for 
tiH production in the large Higgs energy end point is not invariant under separate rotations 
of the spin quantization axis. However, for our purposes it is sufficient to define a triplet 
Wilson coefficient that is averaged over the three spin configurations. Using such an aver- 
aged triplet Wilson coefficient the Higgs energy spectrum at NLL order can also be cast in 
the simple form of Eq. 
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FIG. 2: The singlet form factor Fq and the different contributions to the triplet form factor F^ ,Z for 
run = 120 GeV (respective upper lines), ran = 140 GeV (respective middles lines), and mjj = 160 
GeV (respective lower lines). Here, F^ ,Z 2 refers to the contribution to F^' Z which is proportional 

to Af, and F\'\ tgz an d ^i'^ are defined analogously. The top mass is set to mt = 180 GeV and 
the other parameters are given in Eqs. (|35|). 



H 




FIG. 3: Diagrams describing virtual one-loop QCD corrections in the Standard Model to com- 
pute the NLL matching conditions for the operators Ji jP and Jo iP that describe ti production for 
invariant masses Q 2 ~ Am 2 . Self energy diagrams are implied. 



IV. MATCHING CONDITIONS AT NLL ORDER 

At NLL order, we need to account for the 0(a 2 ) contributions to the Coulomb potential 
in Eq. (jHJ), the NLL running of the coefficients C\ and Co and their 0(a s ) matching conditions 
at v = 1. The latter QCD corrections are process specific and cannot be inferred from results 
obtained in earlier computations for the ti threshold in e + e~ ti. Since we are interested 
in the region where Q 2 m 4m 2 and we can neglect gluon exchange with the Higgs final 
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<5ci 
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175 


120 


—2.266(1) 


—2.339(1) 


—2.397(1) 


—2.453(1) 


—2.509(1) 


—2.564(1) 




140 




-2.339(1) 


-2.401(1) 


-2.459(1) 


-2.516(1) 


-2.572(1) 




160 




-2.340(0) 


-2.407(1) 


-2.467(1) 


-2.524(1) 


-2.581(1) 


180 


120 


-2.254(1) 


-2.329(1) 


-2.386(1) 


-2.441(0) 


-2.495(1) 


-2.548(1) 




140 




-2.327(1) 


-2.389(1) 


-2.446(1) 


-2.501(1) 


-2.555(1) 




160 




-2.328(1) 


-2.394(1) 


-2.452(0) 


-2.509(0) 


-2.564(1) 



TABLE I: Numerical values for the averaged spin triplet matching condition 5ci(\/s, mj, m#) for 
typical values of -y/i, m t and ma- The masses and energies are given in units of GeV. 



5 co 




m t 


m H 


500 


600 


700 


800 


900 


1000 


175 


120 


-0.562(5) 


-0.638(4) 


-0.692(4) 


-0.733(1) 


-0.763(2) 


-0.787(1) 




140 




-0.636(3) 


-0.691(4) 


-0.732(3) 


-0.762(2) 


-0.786(1) 




160 




-0.634(3) 


-0.690(4) 


-0.731(4) 


-0.762(1) 


-0.785(2) 


180 


120 


-0.553(4) 


-0.627(3) 


-0.684(2) 


-0.725(1) 


-0.756(1) 


-0.780(2) 




140 




-0.625(4) 


-0.682(4) 


-0.725(4) 


-0.755(2) 


-0.780(2) 




160 




-0.622(4) 


-0.681(4) 


-0.723(1) 


-0.755(1) 


-0.779(2) 



TABLE II: Numerical values for the spin singlet matching condition 5co(^,mt,mH) for typical 
values of y^, mt and mjf. The masses and energies are given in units of GeV. 

state, we can apply the same nonrelativistic powercounting as the one known for the process 
e + e~ — > ti. Thus real gluon radiation, which is related to emission of gluons with energies of 
order m t v 2 ~ m t a 2 s in the ti cm. frame, is suppressed by a factor a s x {Q 2 — Am 2 ) /m 2 ~ v 3 
(using the nonrelativistic counting v ~ a s ). Therefore we only need to take into account 
the hard contributions contained in the virtual QCD corrections shown in Fig. El We have 
extracted the hard QCD corrections from the results for the Standard Model amplitude for 



e + e~ — > tiH computed in Ref. jljj. The results of Ref. [ll[ were given in term of numerical 
codes for form factors to a number of standard matrix elements. We have computed the 
contributions of these form factors to da/dEn in the large Higgs energy region. To the result 
we have matched the corresponding expression for the 0(a s ) corrections contained in the 
effective theory expression in Eq. PK|) at jj, — mt {p — 1). This allowed us to determine the 
0(a s ) matching conditions to co,i numerically, 

c 0> i(p = 1) = 1 + CFOls ^ m ^ Sc ,i{y/s, m t , m H ) ■ (32) 
With this method the averaged triplet coefficient is obtained automatically. Since the 



numerical results of Ref. are functions of the four-momenta of the initial and final state 
particles we were able to approach the large Higgs energy end point from several directions 
in the final state phase space in order to test the numerical stability and the independence of 
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the result on the direction with which the end point is approached. For example, depending 
on the top-antitop quark cm. relative velocity v, the Higgs and the top and antitop quark 
four momenta (moving in the x-y plane to be definite) can be written as 



p H = (E H , 0, k H , 0) , 
k\ = (Ei, — |ki| sin#, |ki| cos#, 0) , 
k 2 = (E 2 ,-\k 2 \sm(d + 0),\k 2 \cos(d + 0),0) 



(33) 



where 



E; 



m 2 , 



E 



H 



2^~s 



(s + m 2 H - Q 2 ) , Q 2 



4m 2 
1-v 2 



Ex = y/s — Eh — E> 



2 • 



Eo 



s-E H -^ 



x_ (A(E 2 H 



m 



Hi 



(s - 2 y/s E H + m 2 H 
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(34) 



The variable x can be chosen between and 1 and interpolates between the two extreme 
situations x — 0, which refers to the symmetric configuration where Ex = E 2 , and x = 
1, which describes the situation where the top and antitop quarks move into the same 
direction having the maximal possible energy difference. We have extracted the functions 
Scq } x numerically in the limit v — > and estimated the uncertainties from extrapolating to 
v — and from varying the variable x between 0.1 and 0.9. For x too close to zero or one 
the extraction procedure becomes unreliable due to numerical instabilities that develop in 
the code of Ref. 
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in these particular edges of the phase space. In order to identify the 
singlet and the triplet contributions in the full theory we used the helicity basis for the top 
and antitop spinors in the large Higgs energy end point where k x = k 2 to construct the 
singlet and the three triplet tt spin configurations from the standard matrix elements given 
in Ref. 
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The results (including our error estimate) for some representative values for 
y/s, m t and m# are shown in Tab. U for 8cx and in Tab. ITTlfor 5cq. In particular, for smaller 
cm. energies we find that the absolute uncertainties for the singlet matching conditions are 
a few times larger than for the triplet matching conditions. This is a numerical effect that 
can be understood from the fact that for smaller cm. energies the contributions coming from 
the singlet tt spin configuration are much smaller than those coming from the triplet spin 
configuration. In both cases, however, the relative uncertainties are well below 1%. 



V. NUMERICAL ANALYSIS 



In Fig.0]the Higgs energy spectrum in the large energy end point region is displayed at LL 
(dashed lines) and NLL (solid lines) order in the nonrelativistic expansion for the vNRQCD 
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FIG. 4: Higgs energy spectrum in the end point region at LL (dotted lines) and NLL (solid lines) 
order in the nonrelativistic expansion for the vNRQCD renormalization parameters v = 0.1, 0.2, 0.4 
and for ,/i = 600 GeV, mf = 180 GeV, m H = 140 GeV. 



renormalization parameters v = 0.1,0.2,0.4 and for = 600 GeV, m* s = 180 GeV, 
m H = 140 GeV, and 

T t = 1.55 GeV, 

M z = 91.1876 GeV , M w = 80.423 GeV, (35) 
a- 1 = 137.034 , c w = M w /M z . 

At LL order the upper (lower) curve corresponds to v — 0.1 (0.4), while at NLL order the 
upper (lower) curve corresponds to v — 0.2 (0.1). The curves in Fig. 0] show the typical 
behavior of the nonrelativistic expansion for any choice of parameters. While the LL predic- 
tions have a quite large renormalization parameter dependence at the level of several tens of 
percent, the NLL results are stable. Here, the renormalization parameter variation is around 
5%. The stabilization with respect to renormalization parameter variations at NLL order 
arises mainly from the inclusion of the 0(a s ) QCD corrections to the Coulomb potential, 
Eq. (jHJ). Moreover, the NLL curves lie considerably lower than the LL ones. This behavior is 
well known from the predictions for e + e~ — > tt at threshold 3 3 an d originates from the 



structure of the large negative 0(a s ) QCD corrections to the Coulomb potential, Eq. 
and from the sizable negative QCD corrections to the matching conditions in Eq. (J32j) . 

In principle this behavior is a point of concern because it could indicate that the renor- 
malization parameter variation might be an inadequate method to estimate theoretical un- 
certainties. Fortunately, the top quark mass is sufficiently large such that the regions where 
the conventional fixed-order expansion (in powers of the strong coupling) and where the 
nonrelativistic expansion (described by the effective theory) can be applied are expected to 
overlap. In Figs. 03 the Higgs energy spectrum is displayed in the fixed-order expansion at 
the Born (dashed lines) and the 0(a s ) level (solid lines) and in the nonrelativistic expansion 
at LL (dotted lines) and NLL order (solid lines). For the four panels we have chosen the 
parameters (a) yfs = 500 GeV, m H = 120 GeV, (b) ^1 = 550 GeV, m H = 140 GeV, (c) 
^ = 600 GeV, m H = 140 GeV, and (d) ^/s = 700 GeV, m H = 140 GeV. For all cases 
we used m\ s = 180 GeV as the top quark mass. The Born curves were obtained from the 
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FIG. 5: Higgs energy spectrum at LL (dotted lines) and NLL (solid lines) order in the nonrela- 
tivistic expansion for the vNRQCD renormalization parameters v = 0.1,0.2,0.4 and at the Born 
level and at 0(a s ) for ll = ^/s, y/s\v\ for the parameters (a) -y/s = 500 GeV, ran = 120 GeV, 
(b) yfi = 550 GeV, m H = 140 GeV, (c) ^fs = 600 GeV, m H = 140 GeV, (d) </s = 700 GeV, 
run = 140 GeV. For all cases that top mass has been set to m^ s = 180 GeV 



analytic results given in Ref. j7( 4 and the Q (a s ) curves were obtained from the numerical 
program developed by Denner et al. in Ref. [ljj. The two 0(a s ) curves correspond to the 
renormalization scales fi = y/s (lower curve) and /i = \\/sv\ (upper curve), where v is the 
tt relative velocity defined in Eq. (|16J). The latter choice for the fixed-order renormaliza- 
tion scale is motivated by the fact that the relative momentum of the top pair is the scale 
that governs the Coulomb singularities contained in the fixed-order expansion close the large 
Higgs energy end point. This choice for the fixed-order renormalization scale is therefore the 
more appropriate one closer to the Higgs energy end point. For the nonrelativistic expansion 
the renormalization parameters v = 0.1,0.2,0.4 were employed, as in Fig. EJ The results 
in Fig. l^-d demonstrate that there is an overlap between the 0(a s ) fixed-order prediction 
and the NLL nonrelativistic one in the region where the tt relative velocity is approximately 
0.2. (The Higgs energy with v = 0.2 is indicated in each panel by the solid vertical line.) 
The overlap improves for increasing cm. energies or decreasing Higgs masses. This indicates 
that in the overlap regions the higher order contributions summed in the nonrelativistic 
prediction and the higher order relativistic corrections contained in the fixed-order result 



4 We note that the Born expression for Higgs energy distribution given in [8| contains typos |33|: In the 
coefficient function G*2 a factor s is missing for the term (—rrijj + s + sxh)(— p 1 + In the coefficient 
functions G5 and Gq an overall factor of minus one is missing. 
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TABLE III: Collection of cross sections and K factors for various cm. energies and Higgs masses 
and top quark mass m| s = 180 GeV. 

are both small. For smaller cm. energies or increasing Higgs masses, on the other hand, 
the NLL nonrelativistic predictions tend to lie slightly above the 0(a s ) fixed-order results 
(for /x = \fsv) illustrating the impact of the higher order corrections to each type of expan- 
sion. The discrepancy, however, remains comparable to the uncertainties estimated from 
the renormalization parameter variation of the NLL nonrelativistic prediction. We therefore 
conclude that the renormalization parameter variation of the NLL order nonrelativistic pre- 
diction should provide a realistic estimate of the theoretical uncertainties in the large Higgs 
energy region. 

Finally, let us discuss the numerical impact of the nonrelativistic contributions in the 
large Higgs energy region on the total cross section. The results in Figs. El show that 
the summation of the Coulomb singularities and the logarithms of the top quark velocity 
lead to an enhancement of the cross section and that the portion of phase space where 
nonrelativistic effects are important increases for decreasing cm. energy or increasing Higgs 
mass. In Tab. lIIII the impact of the summations is analyzed numerically for various choices 
of the cm. energy and the Higgs mass. For all cases the top quark mass mj s = 180 GeV 
was used. 5 The other parameters are fixed as in Figs.E] and discussed earlier. In the table 
(j(Born) refers to the Born cross section and cr(a s ) to the 0(a s ) cross section in fixed-order 



5 For the 0(a s ) fixed-order numbers that were obtained from the numerical code of Ref. 11] for Fig. [S] and 
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perturbation theory using \x = -Js as the renormalization scale, the choice employed in the 
analysis of Ref. [11]. All cross sections are given in fb units. The term cr(NLL) refers to 
the sum of the 0(a s ) fixed-order cross section for v > 0.2 using \i = y/sv and the NLL 
nonrelativistic cross section for \v\ < 0.2 with the renormalization parameter v = 0.2. The 
numbers for a(NLL) represent the currently most complete predictions for the total cross 
section of the process e + e~ — > tiH as far as QCD corrections are concerned. For the results 
for cr(NLL) we have also given our estimate for the theoretical error. For the fixed-order 
contribution (y > 0.2) we have estimated the uncertainty by taking the maximum of the 
shifts obtained from varying \i in the ranges [t/s, 2-^/i], [y/s,y/s/2] and [y/sv,y/s\; for the 
nonrelativistic contribution in the end point we have assumed an uncertainty of 5% for all 
cases. For the errors given in Tab. IIIII both uncertainties were added linearly. 

The results show that the effect of the summations in the large Higgs energy region is 
particularly important for smaller cm. energies and larger Higgs masses, when the portion 
of the phase space where the nonrelativistic expansion has to be applied is large. Here, 
the higher order summations contained in the nonrelativistic expansion can be comparable 
to the already sizable 0(a s ) fixed-order corrections and enhance the cross section further. 
This is advantageous for top Yukawa coupling measurements for the lower cm. energies that 
are accessible in the first phase of the International Linear Collider experiment. For higher 
cm. energies the effect of the nonrelativistic summations of contributions from beyond 0(a s ) 
is less pronounced and decreases to the one-percent level for cm. energies above 700 GeV, see 
column seven of Tab . 11 111 For all cases, except for very large cm. energies around 1000 GeV, 
however, the shift caused by the terms that are summed up in the nonrelativistic expansion 
exceeds the theoretical error. In the last column we have also displayed the ratio of the 
NLL nonrelativistic cross section for |i>| < 0.2 with v = 0.2 and the 0(a s ) fixed-order 
cross section for v < 0.2 with fi = y/s to illustrate the effect of the summation of QCD 
corrections beyond one loop that have to be carried out in the large Higgs energy region. 
Interestingly, the higher order summations lead to correction factors ranging between about 
1.7 and 1.8 that are only very weakly dependent of the cm. energy and the Higgs mass. This 
fact might prove useful for rough implementations of nonrelativistic ti effects in other high 
energy processes. 

VI. CONCLUSION 

In this paper we have computed the Higgs energy distribution for the process e + e~ — > tiH 
in the large Higgs energy end point region. Due to Coulomb singularities ~ (a s /v) n and 
logarithmic divergences ~ (a s \nv) n that arise at higher order of conventional multiloop- 
perturbation theory in the large Higgs energy end point region, v being the top velocity in the 
ti cm. frame, usual fixed-order perturbation theory breaks down. We use the nonrelativistic 
effective theory vNRQCD to sum both types of singular terms and determine the Higgs 
energy distribution in the end point region at next-to-leading logarithmic order. While the 
Coulomb singularities are summed by solving the nonrelativistic equation of motion for the ti 
Green function, the logarithmic terms are summed using the renormalization group equations 



Tab. Hill we used the top pole mass as an input with the approximation mf olc = Wj S . Numerically the 
actual difference is of order a?, see Eq. ©. We have checked that this approximation is justified within 
the theoretical uncertainties shown in the fifth column of Tab. IIIII 
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in the effective theory. We have compared the effective theor y p redictions for the Higgs 
energy spectrum with the 0(a s ) fixed-order results from Ref. [ll| and can show that the 
nonrelativistic and the fixed-order expansion (with the renormalization scale \x set of order 
\J~sv) have a region of overlap where v is around 0.2. For our next-to-leading logarithmic 
order effective theory prediction we estimate a theoretical error of 5%. The summations lead 
to an additional sizable enhancement of the total cross section compared to the positive 
0(a s ) fixed-order corrections and are particularly important for lower cm. energies (or 
larger Higgs masses) where the portion of the phase space where the ti QCD dynamics 
is nonrelativistic is increased. For higher energies the effects decrease and reach the one 
percent level at cm. energies around 700 GeV. We have given updated predictions for the 
total cross section a(e + e~ — > tiH) combining our effective theory next-to-leading logarithmic 
prediction in the large Higgs energy region with the fixed-order 0(a s ) results determined in 
Ref. jnj. 
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